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1.  Introduction 

Recently  we  had  occasion  to  examine  a  number  of  operational  cloud  seeding 
programs  for  possible  effects  of  silver  iodide  ground  seeding  on  the  rainfall  in 
areas  downwind  from  the  nominal  target  area.  This  preliminary  report  gives  the 
results  of  the  analysis  performed  to  date  on  16  operational  programs  in  11  project 
areas  in  eastern  United  States,  involving  a  total  of  62  seeded  months. 

Cloud  seeding  operations  as  objects  for  statistical  evaluation  are  in  contrast 
to  most  of  the  remaining  papers  in  this  symposium  which  are  concerned  with 
scientific  experiments.  However,  we  feel  that  the  results  obtained  as  well  as  some 
of  the  statistical  techniques  used  may  be  of  sufficient  interest  and  importance  to 
be  discussed  along  with  the  numerous  other  problems  in  weather  modification. 

2.  Analysis  and  results 

A  common  procedure  for  evaluating  the  possible  effects  of  cloud  seeding  on 
the  precipitation  in  a  target  area  is  to  compare  the  actual  precipitation  during  a 
seeded  period  with  that  which  would  have  been  expected  in  the  absence  of  seed¬ 
ing.  This  expected  amount  is  determined  from  a  regression  equation  describing 
the  relationship  between  the  average  precipitation  in  the  target  area  and  the 
average  precipitation  in  a  nearby  control  area  which  presumably  is  unaffected 
by  the  seeding  operation.  This  regression  equation  can  be  determined  from  a 
historical  record  during  a  period  when  no  cloud  seeding  operations  were  known 
to  have  taken  place.  A  measure  of  the  effectiveness  of  cloud  seeding  is  the  anom¬ 
aly  d,  the  difference  between  the  actual  precipitation  (Y)  in  the  target  area  and 
that  expected  (?)  if  no  seeding  had  been  carried  out,  that  is, 

(1)  d  =  Y  —  f. 

Since  our  concern  was  more  with  the  pattern  of  anomalies  outside  and  beyond 
the  nominal  target  area  than  with  the  average  over  the  target,  we  decided  in 
advance  of  any  data  processing  to  modify  the  usual  procedure  by  determining 
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precipitation  anomalies  d  for  individual  stations  both  inside  and  outside  the 
original  target.  Since  the  control  areas  were  usually  to  the  south  or  west  of  target 
area,  we  decided  to  examine  the  precipitation  at  stations  lying  in  a  semicircle 
located  generally  to  the  east  or  northeast  of  the  center  of  the  target  area.  In  some 
cases  it  was  possible  to  include  in  this  extended  target  area  stations  as  far  as  150 
or  more  miles  from  the  target  center. 

A  number  of  difficulties  arose  at  this  stage  that  produced  a  lack  of  uniformity 
in  the  data  available  for  analysis  of  the  11  project  areas.  The  map  in  figure  1 


Location  of  projects  evaluated. 

showing  the  location  of  these  project  areas  indicates  some  of  these  problems.  It  is 
seen  that  some  of  the  areas  are  near  coast  lines  and  target  areas  extended  down¬ 
wind  would  be  over  the  oceans  where  no  data  are  available.  The  extended  target 
areas  also  differed  widely  in  the  number  and  distribution  of  raingages  and  the 
length  of  historical  records  available  for  climatological  analysis.  The  16  programs 
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also  differed  in  respect  to  their  duration,  some  of  the  seeding  operations  running 
for  periods  as  short  as  one  month  or  less.  It  was  recognized  in  advance  that  the 
noise  level  for  individual  projects  might  be  high  and  some  kind  of  a  composite 
picture  would  be  necessary  for  meaningful  results.  We  felt  that  by  combining  all 
the  projects  there  might  be  sufficient  data  to  indicate  in  a  general  way  how  pre¬ 
cipitation  anomalies  were  related  to  distance  from  the  nominal  target.  An  exam¬ 
ination  of  the  climatic  records  for  continuity  of  data,  completeness,  and  so  forth, 
resulted  in  a  total  of  784  stations  for  analysis.  These  were  distributed  among  the 
various  projects  as  shown  in  the  third  column  of  table  I. 


TABLE  I 

Summary  of  Projects,  Number  of  Months  and  Stations  Used,  for  the 
Analysis  of  Precipitation  in  the  Nominal  and  Extended  Target  Areas 


Project 

Locale 

Number  of  Months 
Seeded  Random 

No.  of 
Stations 

Station  Months 
Seeded  Random 

A 

New  Hampshire 

2 

2 

36 

72 

72 

B 

Massachusetts 

1 

1 

12 

12 

12 

C 

Southeastern  New  York 

1 

1 

39 

39 

39 

D 

Connecticut 

1 

1 

82 

82 

82 

E 

New  Jersey 

3 

3 

84 

252 

252 

F 

Eastern  New  York 

1 

1 

121 

121 

121 

G 

Eastern  Pennsylvania 

2 

2 

104 

208 

208 

H 

Southern  Pennsylvania 

10 

10 

102 

1020 

1020 

I 

South  Carolina 

5 

5 

54 

270 

270 

J 

Mid-Potomac 

24 

3 

62 

1488 

186 

K 

North  Carolina 

12 

12 

88 

1056 

1056 

Total 

62 

41 

784 

4620 

3318 

In  order  to  make  the  analysis  comparable  to  some  others  that  had  been  re¬ 
ported,  the  first  step  was  to  make  a  cube  root  transformation  on  all  the  raw 
precipitation  data,  using  monthly  totals  for  individual  stations  and  the  average 
for  the  stations  in  the  control  areas.  It  should  be  understood  that  all  further  dis¬ 
cussion  is  in  terms  of  data  that  have  undergone  a  cube  root  transformation. 

The  next  step  in  the  statistical  analysis  was  to  derive  for  each  project  a  re¬ 
gression  equation  for  each  station  in  the  target  and  extended  target  areas.  This 
equation  was  of  the  form 

(2)  Y  =  a  +  bX 

where  X  is  the  mean  monthly  precipitation  total  for  stations  in  the  control  area 
and  the  coefficients  a  and  b  were  determined  from  the  historical  data  using  the 
observed  monthly  precipitation  totals.  Figure  2  illustrates  the  distribution  of 
the  stations  in  the  control  and  target  areas  for  Project  A  (New  Hampshire). 
Figure  3  illustrates  a  scatter  diagram  showing  a  typical  regression  line  and  the 
transformed  data  (Bradford,  New  Hampshire). 

In  addition  to  deriving  the  regression  coefficients  for  each  station,  the  correla- 
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tion  coefficient  r  was  computed.  Maps  of  r2  were  then  plotted,  as  shown  in  figure 
4.  In  this  figure  the  value  of  r2  for  each  station  of  Project  A  is  plotted  in  terms  of 
station  coordinates  expressed  in  degrees  of  latitude  and  longitude  from  the  center 


74  73  72  71  70  69 


Figure  2 

Location  of  stations  used  in  control, 
target  and  extended  target  areas  for  Project  A. 

of  the  nominal  target  area.  This  procedure  was  followed  for  all  the  projects  to 
facilitate  the  comparison  and  combining  of  results.  The  maps  similar  to  those 
shown  in  figure  4  were  helpful  in  detecting  possible  errors  and  in  addition  have 
an  interest  from  a  meteorological  or  climatological  viewpoint. 
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After  the  regression  equations  were  computed  it  was  possible  to  determine  for 
each  station  and  for  each  seeded  month  an  anomaly 

(3)  d  =  Y  -  f , 

where  Y  is  the  observed  rainfall  index  and  ¥  is  the  rainfall  estimated  from  the 
regression  equation.  Also  determined  was  the  statistic 


where  E  is  the  error  estimate  of  d  as  given  by  classical  regression  theory.  These 


MONTHLY  MEAN  CONTROL  AREA  PRECIPITATION  (CUBE  ROOT) 

Figure  3 

Regression  line  and  data  for  typical  station 
in  Project  A  (Bradford,  N.  H). 

values  were  computed  for  every  station  in  each  project  for  each  seeded  month. 
Thus,  the  basic  experimental  unit  used  here  was  a  project  month;  in  those  cases 
when  the  seeding  operation  ran  for  only  a  fraction  of  a  month  an  adjustment  was 
made  to  make  it  comparable  to  an  entire  month. 

Our  primary  objective  was  to  determine  how  the  anomalies  d  varied  with 
distance  from  the  target  center.  This  was  done  by  classifying  all  the  stations  in  a 
project  according  to  distance,  using  intervals  of  0.4°.  This  class  interval  was  a 
compromise  since  a  smaller  interval  left  too  many  classes  with  small  numbers  of 
stations  and  use  of  a  larger  class  interval  presented  the  danger  of  smoothing 
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■  Figure"4 

Plot  of  r2  for  stations  in  Project  A. 

out  details  that  might  be  of  interest  or  importance.  These  various  distance  classes 
were  called  “rings”  and  defined  as  follows: 

Ring  Distance  in  Degrees 

1  0.000  to  0.399 

2  0.400  to  0.799 

3  0.800  to  1.199 

4  1.200  to  1.599 

5  1.600  to  1.999 

6  2.000  to  2.399 

7  2.400  to  2.799 

8  2.800  to  3.199 

9  3.200  to  3.599 

10  3.600  to  3.999 

11  4.000  and  greater. 

The  results  of  averaging  the  d  anomalies  by  ring,  denoted  here  as  D,  are  shown 
in  table  II  according  to  project  months.  The  last  column  in  the  table  shows  the 
average  anomaly  D  for  the  project  month,  giving  equal  weight  to  each  ring.  The 
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TABLE  II 

Average  Anomalies  D  According  to  Distance  Rings  for  the 
62  Seeded  Project  Months 


Proj¬ 

ect 

Month 

Year 

1 

2 

3 

4 

Rings 

5  6 

7 

8 

9 

10 

11 

Mean 

D 

A 

Nov. 

1964 

.05 

.10 

.20 

.12 

.15 

.09 

.16 

.12 

A 

Nov. 

1957 

.13 

.12 

.14 

.11 

.15 

.14 

.17 

.14 

B 

Nov. 

1964 

.04 

.02 

.06 

.04 

C 

Aug. 

1964 

.08 

.08 

.09 

.09 

.09 

D 

Aug. 

1964 

.04 

.05 

-.04 

-.03 

-.04 

-.01 

-.01 

E 

Oct. 

1964 

.09 

.08 

.12 

.14 

.12 

.07 

.09 

.10 

E 

Nov. 

1964 

-.12 

-.12 

-.09 

-.07 

-.01 

-.06 

-.04 

-.07 

E 

Dec. 

1964 

.05 

.03 

.05 

.07 

.06 

.05 

.14 

.06 

F 

Jul. 

1964 

-.09 

.04 

.18 

.19 

.13 

.03 

.19 

.27 

.28 

.14 

G 

Aug. 

1954 

.12 

.12 

.13 

.21 

.24 

.33 

.19 

G 

Jul. 

1955 

-.00 

-.06 

.01 

.19 

.38 

.58 

.18 

H 

Sep. 

1957 

-.21 

-.27 

-.02 

-.03 

.00 

-.13 

-.04 

-.10 

H 

Oct. 

1957 

-.04 

-.02 

-.02 

-.11 

-.13 

-.11 

-.08 

-.07 

H 

Nov. 

1957 

-.00 

-.02 

-.03 

-.01 

-.05 

-.06 

-.10 

-.04 

H 

Dec. 

1957 

-.01 

.01 

.00 

.02 

.03 

.05 

.07 

.02 

H 

Jan. 

1958 

.08 

.08 

.07 

-.00 

.01 

.03 

.02 

.04 

H 

Oct. 

1954 

.10 

.08 

.00 

-.03 

-.14 

-.19 

-.30 

-.07 

H 

Nov. 

1954 

.01 

.08 

.11 

.13 

.13 

.33 

.38 

.17 

H 

Oct. 

1963 

-.29 

-.23 

-.12 

-.24 

-.19 

-.25 

-.18 

-.21 

H 

Nov. 

1963 

.05 

.01 

.03 

.00 

.02 

.01 

-.01 

.02 

H 

Dec. 

1963 

-.01 

-.02 

-.04 

-.06 

-.10 

-.06 

-.09 

-.05 

I 

Feb. 

1957 

.14 

.14 

.10 

.13 

.08 

.28 

.36 

.33 

.33 

.37 

.35 

.24 

I 

Mar. 

1957 

.08 

.05 

.04 

.11 

.02 

-.07 

-.01 

.09 

.09 

-.01 

.05 

.04 

I 

Apr. 

1955 

-.03 

-.09 

-.05 

.09 

.25 

-.00 

-.14 

-.49 

-.57 

-.32 

-.36 

-.16 

I 

May 

1955 

-.02 

-.01 

.04 

.02 

-.07 

.13 

-.05 

-.04 

.10 

-.08 

-.05 

-.00 

I 

Jun. 

1955 

.38 

.24 

.26 

.10 

-.03 

-.00 

.09 

.09 

.22 

.29 

.24 

.17 

J 

Jun. 

1957 

-.02 

-.02 

-.05 

-.10 

-.09 

-.10 

-.03 

-.13 

-.41 

-.11 

J 

Jul. 

1957 

-.06 

-.04 

-.10 

-.07 

-.11 

-.21 

-.14 

-.14 

-.32 

-.13 

J 

Aug. 

1957 

.02 

.06 

.31 

.39 

.32 

.28 

.15 

.12 

.11 

.20 

J 

Jun. 

1958 

.10 

.15 

.24 

.11 

.03 

.12 

.08 

-.06 

-.12 

.07 

J 

Jul. 

1958 

-.05 

-.04 

-.09 

.01 

.00 

.15 

.10 

-.02 

-.01 

.01 

J 

Aug. 

1958 

.01 

-.01 

.02 

.12 

.14 

.23 

.09 

-.00 

.00 

.07 

J 

Jun. 

1959 

.21 

.24 

.23 

.14 

.23 

.11 

.13 

.08 

.30 

.19 

J 

Jul. 

1959 

.26 

.04 

.09 

.16 

.19 

.20 

.10 

.22 

.15 

.16 

J 

Aug. 

1959 

.04 

.10 

.12 

.14 

.14 

.15 

.20 

.17 

.35 

.16 

J 

Jun. 

1960 

.10 

.23 

.05 

.01 

.04 

-.05 

.00 

.02 

-.21 

.02 

J 

Jul. 

1960 

-.07 

-.01 

.08 

.22 

.29 

.36 

.30 

.24 

.34 

.19 

J 

Aug. 

1960 

.26 

.26 

.31 

.34 

.16 

.12 

.13 

.17 

.16 

.21 

J 

Jun. 

1961 

.10 

.06 

.06 

-.06 

.09 

.14 

.10 

-.05 

-.11 

.04 

J 

Jul. 

1961 

.10 

-.06 

-.01 

-.07 

.21 

.06 

.18 

.23 

.46 

.12 

J 

Aug. 

1961 

-.34 

-.18 

-.10 

.01 

-.05 

-.04 

.00 

.03 

.08 

-.07 

J 

Jun. 

1962 

-.01 

.06 

.03 

.03 

.16 

.17 

.06 

-.04 

-.08 

.04 

J 

Jul. 

1962 

.05 

.01 

-.05 

-.12 

-.25 

-.15 

-.19 

-.16 

-.23 

-.12 

J 

Aug. 

1962 

-.27 

-.07 

-.22 

-.20 

.10 

.20 

.13 

.28 

.27 

.02 

J 

Jun. 

1963 

.21 

.24 

.21 

.18 

.19 

.14 

-.10 

-.16 

-.17 

.08 

J 

Jul. 

1963 

-.08 

-.18 

-.16 

-.30 

-.18 

-.18 

-.05 

-.02 

.15 

-.11 

J 

Aug. 

1963 

-.03 

.03 

.06 

.08 

-.03 

.08 

-.02 

.00 

-.10 

.01 

J 

Jun. 

1964 

-.16 

-.14 

-.32 

-.28 

-.16 

-.20 

-.22 

-.18 

-.16 

-.20 

J 

Jul. 

1964 

.03 

.09 

.02 

-.19 

.10 

.18 

.02 

-.10 

.33 

.05 

J 

Aug. 

1964 

-.10 

-.09 

-.01 

-.13 

-.17 

-.17 

-.16 

-.31 

-.60 

-.19 

K 

Oct. 

1954 

.00 

.09 

.21 

.19 

.47 

.34 

-.01 

.18 

K 

Nov. 

1954 

.24 

.14 

.10 

.04 

.04 

-.00 

-.18 

.05 

K 

Dec. 

1954 

.04 

-.02 

-.01 

-.00 

-.01 

.04 

.06 

.01 

K 

Jan. 

1955 

-.00 

.02 

-.02 

.07 

.09 

.04 

.01 

.03 

K 

Feb. 

1955 

.14 

.17 

.11 

.10 

.12 

.01 

-.09 

.08 

K 

Mar. 

1955 

-.13 

-.09 

.14 

.02 

-.01 

.32 

.43 

.10 

K 

Apr. 

1955 

-.06 

-.13 

-.05 

-.03 

-.07 

.18 

.26 

.01 

K 

May 

1955 

.25 

.13 

.04 

.08 

.03 

-.07 

-.24 

.03 

K 

Jun. 

1955 

.10 

.01 

-.00 

-.08 

-.08 

-.10 

.09 

-.01 

K 

Jul. 

1955 

.11 

.09 

-.19 

-.09 

-.16 

-.16 

-.11 

-.07 

K 

Aug. 

1955 

-.29 

-.07 

-.01 

-.03 

-.07 

-.05 

-.01 

-.08 

K 

Sep. 

1955 

.06 

.09 

.06 

-.01 

.15 

.04 

-.03 

.05 

Mean  D 

.02 

.03 

.04 

.03 

.05 

.06 

.03 

.01 

.02 

.05 

.05 
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TABLE  III 

Number  of  Stations  Used  in  Precipitation  Analysis  for  Each  Project 
According  to  Distance  Class 


Project 

1 

2 

3 

4 

Ring 

5 

6 

7 

8 

9 

10 

11 

Total 

A 

6 

14 

2 

3 

6 

4 

1 

36 

B 

4 

4 

4 

12 

C 

8 

9 

14 

8 

39 

D 

8 

16 

20 

20 

13 

5 

82 

E 

11 

10 

14 

14 

16 

13 

6 

84 

F 

5 

7 

14 

18 

23 

15 

14 

20 

5 

121 

G 

8 

8 

33 

37 

15 

3 

104 

H 

5 

12 

17 

26 

21 

14 

7 

102 

I 

2 

7 

7 

8 

3 

5 

9 

2 

1 

3 

7 

54 

J 

3 

8 

7 

6 

8 

6 

13 

9 

2 

62 

K 

8 

19 

13 

17 

20 

9 

2 

88 

Total 

68 

114 

145 

157 

125 

74 

52 

31 

8 

3 

7 

784 

last  line  in  the  table  shows  the  average  for  each  column  or  ring.  Table  III  shows 
the  number  of  stations  in  each  ring  and  it  is  noted  that  the  projects  differ  widely 
in  terms  of  the  number  of  stations  analyzed  and  the  distance  of  these  stations 
from  the  target  center. 

.12 

.08 

.04 


D  0 
-.04 
-.08 
-.12 

1.0  2.0  3.0 

DIST.  FROM  APPROX.  CENTER  OF  TARGET  (DEGREES) 

Figure  5 

Average  anomaly  D  for  62  seeded  months 
according  to  distance  from  center  of  nominal  target. 
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The  averages  at  the  bottom  of  table  II  are  shown  plotted  in  figure  5  along  with 
the  90  and  95  per  cent  confidence  limits  as  determined  from  the  variances  com¬ 
puted  from  the  project  month  values  given  in  table  II.  Figure  6  shows  the  average 


DIST.  FROM  APPROX.  CENTER  OF  TARGET  (DEGREES) 

Figure  6 

Average  anomaly  D  according  to  distance  from  center 
of  the  nominal  target  and  winter  and  summer  seasons. 

according  to  the  winter  and  summer  seasons.  There  are  indications  of  positive 
anomalies,  some  of  them  statistically  significant,  extending  to  distances  of  ap¬ 
proximately  150  miles. 

One  might  suspect  a  hidden  source  of  bias  or  even  computational  or  program¬ 
ming  errors.  To  give  some  protection  against  these  possibilities,  the  original  plan 
of  the  investigation  included  a  set  of  random  or  fictitious  months  chosen  for  each 
project  for  use  as  a  statistical  control.  Forty  one  such  unseeded  months  were 
chosen  from  periods  in  which  no  seeding  took  place,  and  independent  of  the 
historical  data  from  which  the  regressions  were  derived.  The  results  for  these 
random  or  unseeded  months,  analyzed  in  the  same  way  as  the  operational  seeded 
months,  are  presented  in  table  IV.  Figure  7  shows  the  averages  according  to  dis¬ 
tance  from  the  center  of  the  nominal  target  together  with  the  90  and  95  per  cent 
confidence  limits.  Except  for  Ring  1,  which  has  fewer  stations  than  Rings  2-5, 
there  is  no  suggestion  of  significant  departure  from  the  expected.  In  any  case, 
the  suggestion  is  that  any  bias  in  the  analysis  is  on  the  negative  side,  which 
would  make  the  detection  of  any  real  positive  effects  of  seeding  more  difficult. 
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TABLE  IV 

Average  Anomalies  D  According  to  Distance  Rings  for  the 
41  Unseeded  Project  Months  Analyzed 


Proj 

ect 

Month 

Year 

1 

2 

3 

4 

Rings 

5 

6 

7 

8 

9 

10 

11 

Mean 

D 

A 

Nov. 

1963 

-.05 

-.02 

+  .03 

+  .08 

+.06 

+  .05 

+  .03 

+  .03 

A 

Nov. 

1962 

+.03 

+.09 

+.11 

+  .10 

+.13 

+  .16 

+  .19 

+.12 

B 

Nov. 

1963 

-.03 

-.00 

+.06 

+.01 

C 

Aug. 

1963 

+.07 

+.01 

+.05 

+.12 

+.06 

D 

Aug. 

1963 

-.02 

+.08 

+  .05 

+  .08 

+.09 

+.16 

+.07 

E 

Oct. 

1963 

-.13 

-.11 

-.08 

-.02 

-.02 

-.05 

+  .20 

-.03 

E 

Nov. 

1963 

+.11 

+.08 

+.06 

+.03 

+.02 

+.00 

+  .11 

+.06 

E 

Dec. 

1963 

-.02 

+.00 

-.00 

+  .03 

-.00 

-.03 

-.04 

-.01 

F 

Jul. 

1963 

-.15 

-.15 

+.00 

+.08 

-.01 

+.05 

+.17 

+.15 

-.05 

+  .01 

G 

Aug. 

1956 

-.13 

-.09 

-.11 

-.13 

-.14 

-.13 

-.12 

G 

Jul. 

1957 

-.12 

-.26 

-.23 

-.23 

-.05 

-.00 

-.15 

H 

Jan. 

1960 

-.21 

-.05 

+.07 

-.10 

-.10 

-.09 

+.02 

-.07 

H 

Oct. 

1960 

+.00 

-.02 

+  .04 

+  .07 

+.05 

+.08 

+.05 

+.04 

H 

Nov. 

1960 

-.06 

-.11 

-.09 

-.01 

-.02 

-.10 

-.08 

-.07 

H 

Dec. 

1960 

-.00 

-.02 

-.02 

-.07 

-.07 

-.14 

-.17 

-.07 

H 

Sep. 

1959 

+  .09 

-.02 

-.05 

-.04 

-.07 

-.08 

-.15 

-.05 

H 

Oct. 

1962 

+.26 

+  .14 

+  .06 

-.05 

-.02 

+.05 

+.13 

+.08 

H 

Nov. 

1962 

+.00 

-.08 

-.08 

-.12 

-.19 

-.15 

-.23 

-.12 

H 

Oct. 

1964 

+.10 

+  .05 

-.03 

-.09 

-.12 

-.20 

-.29 

-.08 

H 

Nov. 

1964 

-.10 

-.06 

-.07 

-.10 

-.11 

-.12 

-.16 

-.10 

H 

Dec. 

1964 

+.01 

-.00 

-.01 

+.01 

-.00 

+.02 

-.01 

+.00 

I 

Feb. 

1956 

-.05 

-.06 

-.04 

+.08 

+  .14 

+.09 

+.10 

+.03 

+.20 

+.05 

+.18 

+  .07 

I 

Mar. 

1956 

-.09 

-.02 

-.10 

-.09 

+.02 

-.00 

-.09 

-.01 

-.20 

+.04 

+.03 

-.05 

I 

Apr. 

1956 

-.26 

-.21 

-.09 

-.05 

-.07 

-.04 

-.07 

-.10 

-.18 

-.06 

+.06 

-.10 

I 

May 

1956 

+.14 

-.01 

+.08 

+  .11 

-.04 

+.24 

+.09 

+.35 

+  .42 

-.08 

-.05 

+.11 

I 

Jun. 

1956 

+.00 

+.07 

-.03 

+.04 

+.02 

+.00 

-.07 

-.16 

+.28 

-.14 

+.10 

+.01 

J 

Jun. 

1956 

-.04 

-.02 

-.04 

-.02 

+.01 

+.07 

+.00 

+.06 

-.02 

+.00 

J 

Jul. 

1956 

+.19 

+.24 

+.17 

+.16 

+.05 

+.17 

+  .14 

+.00 

-.09 

+  .11 

J 

Aug. 

1956 

-.10 

-.15 

-.23 

-.45 

-.23 

-.15 

-.12 

-.12 

-.29 

-.20 

K 

Oct. 

1955 

-.18 

-.04 

-.10 

+.03 

-.05 

-.07 

-.05 

-.07 

K 

Nov. 

1955 

-.11 

-.04 

-.01 

-.03 

-.05 

-.09 

-.33 

-.09 

K 

Dec. 

1955 

-.28 

-.24 

-.15 

-.23 

-.32 

-.15 

-.25 

-.23 

K 

Jan. 

1956 

-.16 

-.10 

-.07 

-.04 

-.12 

-.08 

+  .02 

-.08 

K 

Feb. 

1956 

+.01 

+.02 

+.06 

+.07 

+.03 

+.22 

+.30 

+  .10 

K 

Mar. 

1956 

-.10 

-.05 

+.01 

+.06 

-.02 

+.22 

+.37 

+  .07 

K 

Apr. 

1956 

+.05 

+.07 

+  .09 

+.08 

+.00 

+.18 

+  .12 

+.08 

K 

May 

1956 

-.06 

-.01 

+.02 

-.06 

-.02 

-.06 

-.20 

-.06 

K 

Jun. 

1956 

-.18 

-.13 

-.10 

-.10 

-.05 

+.05 

+.20 

-.04 

K 

Jul. 

1956 

+.11 

+.07 

+.06 

+.08 

+.11 

+  .17 

+  .15 

+  .11 

K 

Aug. 

1956 

-.15 

-.14 

+.00 

-.04 

-.03 

-.08 

+.17 

-.04 

K 

Sep. 

1956 

+.17 

+.24 

+.16 

+.13 

+.23 

+.19 

+.07 

+.17 

Mean  D 

) 

-.04 

-.03 

-.01 

-.02 

-.02 

+.01 

+.01 

+.02 

+  .01 

-.04 

+  .06 

Using  the  data  of  table  II  and  table  IV,  it  is  possible  to  make  some  compara¬ 
tive  tests  between  the  seeded  and  unseeded  months.  Table  V  shows  a  comparison 
made  using  the  last  columns  of  tables  II  and  IV,  where  all  stations  in  the  target 
and  extended  target  areas  are  used  irrespective  of  their  distances  downwind. 
The  average  anomaly  for  the  seeded  cases  is  positive  (0.0334)  and  differs  sig¬ 
nificantly  from  zero.  The  unseeded  months  have  a  slightly  negative  but  not 
significant  departure  from  zero.  The  standard  deviations  are  practically  identical 
(0.1103  compared  with  0.0919).  The  standard  t  test  for  differences  between  two 
means  shows  statistical  significance. 

Table  VI  shows  a  comparison  between  the  seeded  and  unseeded  months  by 
individual  project.  In  only  one  project  of  the  eleven  does  the  average  anomaly 
for  the  unseeded  months  exceed  the  anomaly  for  the  seeded  months.  The  bi- 
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DIST.  FROM  APPROX.  CENTER  OF  TARGET  (DEGREES) 

Figure  7 

Average  anomaly  D  for  41  unseeded  months  according  to 
distances  from  center  of  the  nominal  target. 

nomial  sign  test  gives  a  probability  of  p  =  0.01  for  this.  The  unseeded  months 
had  a  mean  anomaly  of  0.00  and  the  seeded  months  had  a  mean  anomaly  of  0.06. 
A  t  test  on  the  paired  comparisons,  with  only  10  degrees  of  freedom  available, 
showed  statistical  significance  on  the  one  tail  test,  but  not  on  the  two  tail  test. 

TABLE  V 


Comparative  Analysis  of  Mean  Precipitation  Anomaly  D  for  Seeded  and 
Unseeded  Months  for  All  Stations  in  Nominal  and  Extended  Target 


Seeded 

Unseeded 

(random  or  control) 

Comparison 

Test 

Number  of  Project  Months 

62 

41 

(101  DF) 

d  =  difference 

0.0334 

-0.0127 

difference  = 

S  =  standard  deviation 

0.1103 

0.0919 

0.0461 

0.1034 

E  =  standard  error  of  d 

0.0140 

0.0144 

‘-J 

2.39 

-0.88 

2.22 

probability  (two  tails) 

p  =  0.020 

p  =  0.40 

p  =  0.03 

probability  (one  tail 
for  increase) 

p  =  0.010 

p  =  0.80 

p  =  0.015 
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TABLE  VI 

Mean  Precipitation  Anomaly  (Cube  Root  Scale)  by  Project  Area 
for  All  Stations  in  Nominal  and  Extended  Target  Areas 


Project 

Locale 

Seeded 

Unseeded 
(random  or 
control) 

Seeded 

Minus 

Unseeded 

A 

New  Hampshire 

0.13 

0.07 

0.06 

B 

Massachusetts 

0.04 

0.01 

0.03 

C 

Southeastern  New  York 

0.09 

0.06 

0.03 

D 

Connecticut 

-0.01 

0.07 

-0.08 

E 

New  Jersey 

0.04 

0.01 

0.03 

F 

Eastern  New  York 

0.14 

0.01 

0.13 

G 

Eastern  Pennsylvania 

0.18 

-0.13 

0.31 

H 

Southern  Pennsylvania 

-0.03 

-0.04 

0.01 

I 

South  Carolina 

0.06 

-0.01 

0.07 

J 

Mid-Potomac 

0.03 

-0.03 

0.06 

K 

North  Carolina 

0.02 

-0.01 

0.03 

Means 

0.06 

0.00 

0.06 

tlO  DF  = 

1.99,  two  tail 

p  =  0.074, 

one  tail  p  —  0.037 

The  results  of  the  analysis  to  date  of  the  average  precipitation  anomalies  in  the 
target  and  extended  target  areas  seem  to  indicate  the  following: 

(1)  the  seeded  months  are  different  from  the  historical  months; 

(2)  the  seeded  months  are  different  from  the  unseeded  months  (chosen  inde¬ 
pendent  of  the  historical  record) ; 

(3)  the  independent  unseeded  months  show  slight  but  not  generally  statis¬ 
tically  significant  anomalies  in  the  negative  direction,  suggesting  that  if  any 
bias  exists  in  the  statistical  analysis,  it  tends  to  operate  in  a  direction  against 
finding  positive  effects  of  seeding. 

If  these  results  had  been  obtained  from  scientific  experiments,  one  would  ordi¬ 
narily  conclude  that  cloud  seeding  had  affected  rainfall.  However,  since  these 
were  not  randomized  experiments,  a  rigorous  statement  regarding  a  causal  rela¬ 
tion  between  seeding  and  precipitation  anomalies  cannot  be  made.  There  may 
be  hidden  biases  and  we  are  continuing  to  examine  the  data  with  such  possibilities 
in  mind.  For  example,  we  plan  to  examine  the  terminal  month  of  the  operations 
to  see  whether  there  was  any  tendency  to  stop  the  seeding  when  higher  than 
normal  precipitation  occurred  in  the  area.  Also,  we  expect  to  examine  the  pre¬ 
cipitation  anomalies  during  the  period  after  the  cessation  of  seeding  since  some 
investigations  have  reported  evidence  of  a  persistence  effect  of  silver  iodide  due 
to  contamination  of  trees,  vegetation,  and  so  on. 

Another  type  of  bias  might  be  introduced  if  cloud  seeding  operations  and/or 
experiments  are  started  in  drought  periods.  This  possibility  has  been  suggested 
by  several  investigators  although  no  quantitative  estimates  have  been  given  re¬ 
garding  the  direction  or  magnitude  of  such  possible  effects.  Our  analysis,  based 
on  the  historical  data  from  unseeded  periods,  disclosed  no  evidence  that  the  re- 
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gression  relationship  between  target  and  control  areas  is  disturbed  during 
drought  periods. 

The  results,  shown  in  table  VII,  indicate  that  the  anomalies  d  expected  after 

TABLE  VII 

Average  Anomalies  D  for  the  First  Three  Distance  Rings  According  to 
Elapsed  Time  from  the  Initial  Month  Chosen  on  the  Basis 
of  Drought  in  the  Target  Area 
*  indicates  beyond  5%  significance  level;  **  beyond  1%  significance  level 


Time 


Initial  Month 

1  Month  Later 

1  Year  Later 

2  Years  Later 

N 

D 

t 

N 

D 

t 

N 

D 

t 

N 

D 

t 

Extreme 

Ring  1 

-0.13 

-7.68** 

29 

-0.00 

-0.15 

45 

0.02 

1.09 

41 

0.01 

0.72 

Drought  in 

Ring  2 

50 

-7.22** 

31 

0.02 

1.03 

45 

0.01 

0.70 

40 

0.03 

1.31 

Target  Area 

Ring  3 

-7.27** 

31 

-0.03 

-1.73 

46 

0.01 

0.38 

40 

0.02 

0.80 

Moderate 

Ring  1 

166 

-8.27** 

111 

0.01 

0.14 

148 

1.02 

134 

0.01 

1.00 

Drought  in 

Ring  2 

166 

-7.35** 

111 

-0.00 

-0.29 

151 

0.21 

135 

-0.00 

-0.04 

Target  Area 

Ring  3 

166 

UM 

-8.11** 

108 

-0.02 

-2.03* 

146 

130 

-0.00 

-0.63 

Mean 

-0.09 

-7.65** 

-0.00 

-0.50 

0.01 

0.45 

0.01 

0.52 

a  dry  period  are  essentially  the  same  as  in  the  entire  sample.  Extreme  drought 
months  were  defined  as  those  falling  in  the  lowest  (driest)  one  tenth  of  the  pre¬ 
cipitation  frequency  distribution.  If  the  month  fell  in  the  lower  one  third  of  the 
distribution,  it  was  classified  as  a  moderate  drought.  The  initial  months  must 
show  significance  since  they  were  selected  as  having  deficient  precipitation. 
If  there  is  any  bias,  it  tends  to  be  in  the  negative  direction,  making  it  more  diffi¬ 
cult  to  detect  positive  seeding  effects.  Another  bias  which  would  tend  to  act  in 
the  same  direction  is  due  to  the  contamination  of  the  control  area  by  the  ground 
seeding.  Such  an  effect  might  be  minimized  by  using  a  different  control  area  for 
each  wind  direction,  storm  type,  season,  and  so  forth.  This  should  also  increase 
the  power  of  a  test  of  seeding  effects.  We  estimate  that  by  using  the  methods 
reported  in  this  analysis,  a  total  of  20  seeded  months  would  be  needed  to  give  a 
90  per  cent  chance  of  detecting  a  13  per  cent  increase  at  the  5  per  cent  significance 
level  for  Project  A  (New  Hampshire,  winter).  There  is  good  reason  to  think  that 
a  considerable  improvement  can  be  made  on  this  by  stratifying  according  to 
relevant  meteorological  parameters  and  by  sharpening  up  the  statistical  esti¬ 
mates  by  multiple  regression  methods,  data  transformations,  and  so  forth.  These 
and  other  problems  are  being  investigated  not  only  in  order  to  obtain  a  better 
understanding  of  what  might  have  happened  during  these  particular  operations 
but  as  part  of  a  general  program  of  searching  for  more  powerful  evaluation 
techniques. 


